NEAR-RESONANT, STEADY MODE INTERACTION: PERIODIC, 
QUASI-PERIODIC AND LOCALIZED PATTERNS* 

MARIA HIGUERAt, HERMANN RIECKE*, AND MARY SILBER§ 

Abstract. 

Motivated by the rich variety of complex periodic and quasi-periodic patterns found in systems 
such as two-frequency forced Faraday waves, we study the interaction of two spatially periodic modes 
' that are nearly resonant. Within the framework of two coupled one-dimensional Ginzburg-Landau 

^ equations we investigate analytically the stability of the periodic solutions to general perturbations, 

including perturbations that do not respect the periodicity of the pattern, and which may lead 
to quasi-periodic solutions. We study the impact of the deviation from exact resonance on the 
destabilizing modes and on the final states. In regimes in which the mode interaction leads to 
traveling waves our numerical simulations reveal localized waves in which the wavenumbers are 
resonant and which drift through a steady background pattern that has an off-resonant wavonumber 
C , _ , J . ratio. 



o 



o 

Oh 

d 



1. Introduction. Pattern- forming instabilities lead to an astonishing variety of 
spatial and spatio-temporal structures, ranging from simple, periodic stripes (rolls) 
■ to spatially localized structures and spatio-temporally chaotic patterns. Even within 

the restricted class of steady, spatially ordered patterns a wide range of patterns 
have been identified and investigated beyond simple square or hexagonal planforms 
including patterns exhibiting multiple length scales: supcrlatticc patterns, in which 
the length scales involved arc rationally related rendering the pattern periodic albeit 
on an unexpectedly large length scale, and quasipatterns, which are characterized by 
incommensurate length scales and which are therefore not periodic in space. These 
I more complex two-dimensional patterns have been observed in particular in the form 

. of Faraday waves on a fluid layer that is vertically shaken with a two-frequency peri- 

' odic acceleration function ^ 1^31 12] , and to some extent also in vertically vibrated 

. Rayleigh-Benard convection 01] and in nonlinear optical systems The Faraday 

system is especially suitable for experimental investigation of pattern formation in 
. systems with two competing spatial modes of instability since such codimension-two 

points are easily accessible by simply adjusting the frequency content of the periodic 
forcing function |16l [S] . The observation of both superlattices and quasipatterns in 
this physical system raises an intriguing question concerning the selection of these 
kinds of patterns: what determines whether, for given physical parameters, a periodic 
or a quasi-periodic pattern is obtained? This provides the main motivation for the 
present paper in which we address certain aspects of the selection problem within the 
^ ' somewhat simple framework of mode competition in one spatial dimension. 

To capture the competition between commensurate/incommensurate length scales 
we focus on the interaction between two modes with a wavelength ratio that is close 
to, but not necessarily equal to, the ratio of two small integers. We are thus led 



(N 
> 



> 

X 



*The work of MH was supported by the Spanish Ministerio de Ciencia y Tecnologia (BFM2001- 
2363). The work of HR was supported by NASA grant NAG3-2113, the Department of Energy 
(DE-FG02-92ER14303) and NSF (DMS-9804673). The work of MS was supported by NASA grant 
NAG3-2364, and by NSF under DMS-9972059 and the MRSEC Program under DMR-0213745. 

tE. T. S. Ingenieros Aeronauticos, Universidad Politecnica de Madrid, Plaza Cardenal Cisneros 
3, 28040 Madrid, Spain. (mariaOfmetsia. upm.es). 

■^Department of Engineering Sciences and Applied Mathematics, Northwestern University, 
Evanston, IL 602028, USA. 

§ Department of Engineering Sciences and Applied Mathematics, Northwestern University, 
Evanston, IL 602028, USA. 



1 



to consider a near-resonant mode interaction. We focus on the case in which both 
modes arise in a steady bifurcation. While at first sight it may seem that an analysis 
of steady-state modes would not be applicable to two- frequency forced Faraday waves, 
for which most of the patterns of interest in the present context have been observed, 
it should be noted that very close to onset the amplitude equations for these waves 
can be reduced to the standing-wave subspace, in which the waves satisfy equations 
that have the same form as those for modes arising from a steady bifurcation (see, for 
example, |35|). 

The competition between commensurate and incommensurate steady structures 
has been addressed previously in the context of a spatially-periodic forcing of pat- 
terns |27[ [71 [TUl I12L I31j . Using a one-dimensional, exiernaZ forcing of two-dimensional 
patterns in electroconvection of nematic liquid crystals, localized domain walls in the 
local phase of the patterns were observed if the forcing wavenumber was sufficiently 
incommensurate with the preferred wavenumber of the pattern j27| . Theoretically, 
the domain walls were described using a one-dimensional Ginzburg-Landau equation 
that included a near-resonant forcing term, which reflects the small mismatch between 
the forcing wavenumber and the wavenumber of the spontaneously forming pattern 
[3 EI. The situation we have in mind in the present paper is similar to these studies 
in so far as the competing modes we investigate also provide a periodic forcing for 
each other. The case of external forcing is recovered if one of the two modes is much 
stronger than the other and consequently the feedback from the forced mode on the 
forcing mode can be ignored. We do not restrict ourselves to this case, however, and 
thus both modes are active degrees of freedom and the interaction between them is 
mutual. 

The analysis of the interaction of two exactly resonating modes near a co-dimension- 
two point at which both modes bifurcate off the basic, unpatterned state has revealed 
a wide variety of patterns and dynamical phenomena. Rich behavior has been found 
in small systems in which the interacting modes are determined by the symmetry and 
shape of the physical domain (e.g. IH EH CH IHIISI EOl ES| ) ■ More relevant for 
our goal are the studies of mode-interaction in the presence of translation symmetry, 
since they allow the extension to systems with large aspect ratio, which are required to 
address the difference between commensurate and incommensurate structures. There- 
fore our investigation builds on a comprehensive analysis, performed by Dangelmayr 
of the interaction between two resonant spatial modes in 0(2)-symmetric systems 
with wavenumbers in the ratio m : n, m < n (Here the 0(2)-symmetry is a conse- 
quence of restricting to spatially periodic patterns in a translation-invariant system.). 
For m > 1 there are two primary bifurcations off the trivial state to pure modes 
followed by secondary bifurcations to mixed modes. These mixed modes can in turn 
undergo a Hopf bifurcation to generate standing wave solutions and a parity-breaking 
bifurcation that produces traveling waves. For m = 1 similar results are found except 
that there is only one pure-mode state, the one with the higher wavenumber; the other 
primary bifurcation leads directly to a branch of mixed modes. Further analyses of 
this system have revealed structurally and asymptotically stable heteroclinic cycles 
near the mode interaction point when m : n — 

1 : 2 El ESI El- More recently, complex 
dynamics organized around a sequence of transitions between distinct heteroclinic cy- 
cles has been discovered in resonances of the form 1 : n I34j . in particular, in the 
cases n = 2 and n = 3. 

Although previous work on resonant mode interactions considered only strictly 
periodic solutions it provided insight into various phenomena that were observed 
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experimentally in large-aspect ratio systems involving large-scale modulations of the 
patterns. For example, in steady Taylor vortex flow it was found experimentally 
that not too far from threshold the band of experimentally accessible wavenumbers is 
substantially reduced compared to the stability limit obtained by the standard analysis 
of side-band instabilities in the weakly nonlinear regime ^J] . The origin of this strong 
deviation was identified to be a saddle-node bifurcation associated with the 1 : 2 
mode interaction PHI ■ In directional solidification jlSl localized drift waves have been 
observed, which arise from the parity-breaking bifurcation [21101^1 that is associated 
with the resonant mode interaction with wavenumber ratio 1 : 2 |26| . Subsequently 
such waves have also been obtained in a variety of other systems including directional 
viscous fingering |S7] , Taylor vortex flow 0^ ^Uj , and premixed flames 0] . In our 
treatment of the near-resonant case the mode amplitudes are allowed to vary slowly 
in space. It therefore naturally incorporates phenomena like the localized drift waves 
and the modification of side-band instabilities by the resonance. 

In this paper we study the interaction of two nearly-resonant modes in a spatially 
extended, driven, dissipative system. Near onset we model the slow dynamics of such 
systems by two amplitude equations of Ginzburg-Landau type, one for each mode. 
We focus on the weak resonances (m -I- n > 5) in order to avoid some of the specific 
features of the strong-resonance cases (for example, the structurally stable heteroclinic 
cycles in the case m : n = 1 : 2). Our primary goal is to investigate the transitions 
between periodic and quasi-periodic states that take place as the result of side-band 
instabilities, with an eye on how the detuning from exact spatial resonance influences 
this process. We find that the detuning can play an important role in the selection 
of the final wavenumbers of the modes involved. For example, it can favor a periodic 
to quasi-periodic transition that would otherwise (i.e., in the case of exact resonance) 
result in a second periodic state. Among the various quasi-periodic states that we 
find in numerical simulations are several that consist of drifting localized structures 
with alternating locked (periodic) and unlocked (quasi-periodic) domains. 

It should be noted that at present there is no rigorous justification for the descrip- 
tion of quasi-periodic patterns using low-order amplitude expansions. In fact, the lack 
of straightforward convergence of such an expansion has been investigated recently 
for two-dimensional quasi-patterns 021 ■ We will not discuss these issues; instead we 
use the coupled Ginzburg-Landau equations as model equations that are known to 
be the appropriate equations for periodic patterns and at the same time also allow 
quasi-periodic solutions. 

The organization of the paper is as follows. In Section 2 we set up the coupled 
Ginzburg-Landau equations that are based on the truncated normal form equations 
for the m : ri-resonance. In Section 3 we utilize and build upon on the detailed 
results of Dangelmayr jT3] to describe the stability properties of steady spatially 
periodic states (i.e, pure and mixed modes) with respect to perturbations that preserve 
the periodicity of the pattern. In Section 4 we turn to the question of stability of 
the steady periodic solutions with respect to side-band instabilities. Here we also 
determine how these instabilities are affected by the detuning from perfect resonance. 
Numerical simulations of the system are carried out in Section 5 to investigate the 
nonlinear evolution of the system subsequent to the side-band instability discussed in 
the previous Section. Our concluding remarks are given in section 6. It should be 
pointed out that a very recent preprint by Dawes et al. presents a complementary 
investigation near the 1 : 2 resonance |14| . 
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2. The amplitude equations. We consider driven, dissipative systems in one 
spatial dimension that are invariant under spatial translations and reflections. We 
further assume that in the system of interest there are two distinct spatial modes that 
destabilize the basic homogeneous state nearly simultaneously. The wavenumbers of 
these two modes, qi and 52, correspond to minima of the neutral curves and are 
assumed to be in approximate spatial resonance, i.e., 

ri(qi +£7) ?«g2, (2.1) 

where m and n > m are positive co-prime integers and the term ej measures the 
deviation from perfect resonance. Physical fields, near onset, may then be expanded 
in terms of these two spatial modes: 

u{x, t) = e[Ai (X, r)e*(«i+^T')^ + A2{X, T)e*«=^ + c.c] + • • • (2.2) 

The two modes are allowed to vary on slow spatial and temporal scales X = ex and 
T ~ e'^t, respectively. Note that in the expansion H2.2() we do not expand about about 
the minima qi^2 of the neutral stability curves, which for 7 7^ are not in spatial 
resonance, but take instead a mode Ai which is in exact spatial resonance with A2. 
This choice of Ai and A2 simplifies the equations, avoiding any explicit dependence 
on the spatial variable in the resulting amplitude equations. Furthermore, we allow 
for a small offset between the critical values of the forcing amplitudes Fic of the two 
modes (see Fig. I2.1|l . 
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Fig. 2.1. Sketch of the neutral stability curve for the two modes. 

The equations governing the evolution of Ai and A2 must be equivariant under 
the symmetry operations generated by spatial translations (T^) and spatial reflections 
(R), which act on (^1,^2) as follows: 

T^:(Ai,^2)^(e^™'^Ai,e™^A2), for ^e[0,2^), (2.3) 
R:{Ai,A2)^ {AuA2), 

where the bar denotes the complex conjugate. Consistent with this equivariance 
requirement the slow evolution of Ai and A2 can be approximated, after rescaling, by 
the Ginzburg-Landau equations 
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A^T = Ml + SA^xx - ilA^x - {-Ml? + p|^2|')Ai + z^^r^A™, (2.4) 

A2T = (a^ + A//)A2 + 5'A2xx - (s'|A2|' + pViP)^2 + (2.5) 

The subscripts indicate partial derivatives with respect to X and T . The main con- 
trol parameter fj, (x {F — i^ic)/e^ measures the magnitude of the overall forcing. 
In addition, we keep track of the offset in the two critical forcing amplitudes with 
A/i = [Fic — F2c)/£^ (see Fig. 12. and capture the detuning between qi and mq2/n 
with 7 = 2Sj. The local curvature of the neutral stability curves near qi and 52 
is measured by S and S' , respectively. We further assume that the nonlinear self- 
and cross- interaction coefficients satisfy the non-degeneracy conditions ss' ^ and 
ss' — pp' ^ and perform a simple rescaling such that s ~ ±1, s' ~ ±1. 

One goal of this paper is to gain insight into the difference between periodic 
and quasi-periodic patterns in systems with two unstable wavenumbers. If the two 
wavenumbers are not rationally related and their irrational ratio is kept fixed as 
the onset for the two modes is approached, then only terms of the form ylij^jp', 
Z = 1, 2, 3 . . ., appear in the equation for A^. For a rational ratio, however, additional 
nonlinear terms arise, which couple the otherwise uncoupled phases of the two modes 
Ai. For the m : n-resonance the leading-order resonance terms are given by A"~^A1^ 
and A'^~^Ai in the equations for Ai and A2, respectively. In order to explore the 
connection between the rational and the irrational case we consider a wavenumber 
ratio which may be irrational, but its deviation from the ratio m : n is of 0{e). This 
allows the mismatch between the two wavenumbers to be captured by the slow spatial 
variable X and to describe periodic and quasi-periodic patterns with the same set of 
equations H2. 412. 5(1 . Equivalently, we could have expanded in the irrationally related 
wavenumbers qi and q2 associated with the minima of the neutral curves in Fig. 12.11 
Then the resonance terms would introduce space-periodic coefficients with a period 
that is related to the mismatch of the wavenumbers. Our choice of the expansion 
wavenumbers (cf. (|2.1() ') removes this space-dependence and introduces the first-order 
derivative —ijdxAi in its place. 

We focus on the weak resonances, m + n > 5, in which the resonant terms are of 
higher order. We neglect, however, non-resonant terms of the form AijAjl^ j = 1, 2 
and 4 < p £ N) which may arise at lower order. This is motivated by the observation 
that such terms do not contribute any qualitatively new effects for small amplitudes. 
The resonant terms, in contrast, remove the unphysical degeneracy that arises when 
the phases are left uncoupled and can therefore influence dynamics in a significant 
way despite appearing at higher order. Note, however, that near onset the resonant 
terms are typically small and the phase coupling between Ai and A2 occurs on a very 
slow time scale. The coupling becomes stronger further above onset where the weakly 
nonlinear analysis may no longer be valid. 

It is often useful to recast Eqs. H2.4I2.5|I in terms of real amplitudes Rj > and 
phases by writing Aj = Rje^"^^ . This leads to 

RiT = pRi - {sRl + pRl)Ri + i^R"^^R^ cos(n^i - m^a) 

+SR1XX - HlxRi - 270ixi?i, (2.6) 
R2T = (m + ^^^)R2 - {s'Rl + p'Rl)R2 + v'R'^-^Rl cos(n0i - m^s) 

+5'R2xx - 5'^lxR2. (2.7) 
RkI^it = -iyRi-^RTsm{n(l)i ~ m(f>2) + HixxRi + 2S(f>ixRix + 2jRix, (2.8) 
R2(f>2T = ly'RT'^Ri sin(n0i - m02) + S'(t)2xxR2 + 2(5>2xi?2X. (2.9) 
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If the spatial dependence in Eqs. H2.6I2.9|I is ignored, the system reduces to a set of 
ordinary differential equations (ODEs) , equivalent to the one analyzed by Dangelmayr 
In this simplified problem, the translational symmetry (T^) causes the overall 
phase to decouple and leaves only the three real variables R2, and the mixed 
phase 

(/) = n0i-m02, (2.10) 

with dynamically important roles. Dangelmayr's bifurcation analysis produced ex- 
pressions for the location of primary bifurcations to pure mode solutions, secondary 
bifurcations to mixed mode solutions, and, in some instances, tertiary bifurcations to 
standing-wave and traveling-wave solutions. These results apply to a general m ; n 
resonance, and prove useful in what follows. 

3. Steady Spatially-Periodic Solutions. In this section we analyze steady 
solutions of Eqs. H2.4I2.5|I of the form 

where i?i_2 > 0, and ^1.2 and k are real. Such states represent spatially periodic 
solutions of the original problem with wavenumbers qi = qi + £7 + ek and 92 = 
q2 + enk/m so that qin = q2m- These solutions break the continuous translational 
symmetry (T^) but remain invariant under discrete translations. 

Within this family of steady states there are generically only two types of non- 
trivial solutions, pure modes and mixed modes, which we describe below. 

I. Pure modes (51^2). These are single-mode states, which take one of two forms: 

Si : {Ai,A2) = {^/^ e'^''^+'^'\0), for m > 1, 
^2 : (Ai, ^2) = (0, e*(("fc/™)^+02))^ (3^2) 

where (j>i,(j>2 £ [0, Stt) and 

a ~ ji — Sk^ + 7fc, 

13 ^ ^1 + 1^^1-5' [nk/mf. (3.3) 

Note that pure modes of type Si are not present if m = 1 (see Eqs. ()2 .412 .5(1 '). More- 
over, the pure modes S'l and S2 are not isolated, but emerge as circles of equivalent 
solutions (parametrized by ^1 or (j)2)- Hereafter we consider resonances m : n where 
m > 2, in which case both pure modes S'l and S2 are present. 

II. Mixed modes (5±). There are two types of mixed modes. 



S± : {Ai,A2) = (i?i e''^''^+^'\R2 e'iir,k/m)x+i^)^^^ ^3 4) 

satisfying 

S± : cos(0) = ±1, 

{s'a - pp) = {ss' - pp')Rl ± [s'vRl - v' pRDK^'' R^^ ,^ 
(s/3 - p'a) = {ss' - pp')Rl ± {sv'Rl - iyp'R^)R^-^R'l'^. ^ ' 

Here 4> is the mixed phase given by Eq. H2.10() . and a and (3 are defined by Eq. H3.3(l . 
As in the case of the pure modes Si^2, translational symmetry implies that there are 
circles of equivalent mixed-modes states (parametrized by ^1, say). Like the pure 
modes, the mixed modes are invariant under reflections (i?) through an appropriate 
origin. 
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3.1. Stability under homogeneous perturbations. The stability of 5*1 2 and 
S± under homogeneous perturbations can be obtained from Dangelmayr's analysis 

, which we review here in some detail to provide the background necessary for our 
analysis. The stability regions in the (a, /9)-unfolding plane simply need to be mapped 
to the (fc, /i)-plane with the (nonlinear) transformation (|3.3() . Each intersection of the 
curves a = and f3 — corresponds to the codimension-two point of |13) : there are 
generically zero or two such intersections. Since the nonlinear coefficients are identi- 
cal in the vicinity of both intersections, the response of the system to homogeneous 
perturbations in corresponding neighborhoods of the (two) intersections is identical. 
In particular, any bifurcation set arising from one intersection arises from the other 
as well. Several examples are given in Fig. 13.11 which shows the various bifurcation 
sets in the (fc, ^)-plane in four representative cases. These bifurcations are described 
below. Note that for mode Ai the deviation of the wavcnumber from the critical 
wavenumber is given by sk, whereas for mode A2 it is given by enk/m. Additional 
details, valid in sufficiently small neighborhoods of the intersections, are available in 

In this paper we study H2.4I2.5|I with the coefficients u and i^' of the higher- 
order resonance terms taken to be of order 1. The resulting stability and bifurcation 
diagrams therefore contain certain features that do not remain local to the bifurcation 
when v' — > 0, i.e. in this limit these features disappear at infinity, fi —>■ 00, and do 
not represent robust aspects of the mode-interaction problem. We return to this issue 
briefly at the end of this section, and indicate which features of our sample bifurcation 
sets are not robust in the limit v,!/' ^ 0. 

The pure modes Si and 6*2 bifurcate from the trivial state when a = and 
(3 = 0, respectively: the bifurcation to Si {S2) is supercritical if s — 1 (s' = 1). 
Their stability is determined by four eigenvalues, one of which is forced to be zero by 
translation symmetry. In the case of Si the remaining three eigenvalues are 

a'^^ = -2a/s, 

Li : X'i'> = (/3 ~ p'a)/s ± \v' , (3.6) 

where the bracketed term with subscript m = 2 is present only if m = 2. For S2 the 
eigenvalues are given by 

L2: \f =\^^^ = {s'a-p(i)/s'. (3.7) 

When one of the eigenvalues {^^±) changes sign the pure modes Si (S'2) 
become unstable to mixed modes, respectively. For S'2 the two eigenvalues coincide 
and the bifurcation occurs along a single line denoted by L2 in Fig. 13. II Similarly, for 
m > 2 the transition from Si to the mixed modes 5*^ occurs along the single curve Li 
fFig. l3.lb ,.b'l. For m = 2 the eigenvalues A^'' are not degenerate and the line Li splits 
into two curves, L'^ and L^; the mixed mode 5+ bifurcates at and S- bifurcates 
at [Lf in Figs. l3.lb -d'). The splitting is due to the presence of the resonant terms 
which are linear in yl2 if m = 2 but not otherwise. Both curves, and , become 
tangent to each other at the intersection a = (3 = 0. 

The response of the mixed modes to amplitude perturbations is decoupled (due to 
reflection symmetry) from the effect of phase perturbations. The amplitude stability 
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Fig. 3.1. Bifurcation sets indicating the creation of pure modes (thick solid lines), mixed modes 
(thin solid lines , L\ and L2), standing waves (dashed lines SW^) and traveling waves (dash- 
dotted lines TW^), as well as saddle-node bifurcations (dashed lines SN). The specific resonance 
m : n is indicated above the plots. In each case Afi = 0.5, 7 = 0.5, S = 5' = I, p = 0.4, and 
p' = 0.67. In (a) s = -s' = 1, u = 0.62 and v' = 1.02, but in (b), (c) and (d) s = s' = -I, v = 0.62 
and v' = -1.02. 



is determined by the eigenvalues of a 2 x 2-matrix M± , whose determinant and trace 
can be written as (recall we consider m > 1) 

det(M±) = -4(p/9' - ss')R\rI ± R'^-^R^'-^H{Ri, R2), (3.8) 
Tr(M±) = -2s(i?? + ss'Rl) ± R^^-^R'^-^ii^in - 2)RI + v'{m - 2)^^). (3.9) 

with 

H{Ri,R2) = -2 [iy's{m - 2)R\ + vs'{n - 2)R^ - RfRKp'iym + pv'n)] (3.10) 

Here Ri and R2 are solutions of Eqs. (|3.5(l . Since in general v and v' are of C'(e™+"~''), 
the contribution H{Ri,R2) from the resonance term affects the steady bifurcation 
determined by H3.8|l only for m < 3. In particular, in the cases m : n ^ 3 : n 
and m : ri = 2 : 3 the function H can balance the first term along curves through 
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the codimcnsion-2 point along which i?2 «C Ri and i?i <i; i?2, respectively, and 
one of the mixed states experiences a saddle-node bifurcation (curves labeled SN 
in Figs. I3.lb .dl \1'6\ . For 2 : n resonances with n > 5 the term Rf drops out in 
H{Ri,R2). Consequently, the i/-term cannot balance the first term in H3.8|l and no 
saddle-node bifurcations occur. 

For m > 4 the sign of det(M±) does not depend on i?i.2. Then the S± solutions 
are always unstable for {pp' — ss') > 0, while for {pp' — ss') < their stability must be 
deduced from the sign of Tr(M±). If ss' = 1 it follows from Eq. H3.9|) that for e ^ 1 
sign(Tr(M-|-)) = — s; the mixed modes S± are then stable to amplitude perturbations 
if s = -1-1. On the other hand, if ss' = —1, the trace of M± can change sign, indicating 
the possibility of a Hopf bifurcation. The resulting time-periodic solutions inherit the 
reflection symmetry of S± (so (t>i ~ 4>2 = 0) and therefore correspond to standing 
waves {SW curves in Fig. 13. lb ). 

Instabilities associated with perturbations of the mixed phase (I2.10|) lead to bi- 
furcations breaking the reflection symmetry. The relevant eigenvalues are given by 

TW : e±^T{i^nR:^ + iy'mRl)R'l-^R^-^, (3.11) 

and may pass through zero only if sign(i^i'') — —1. In this case the S± states undergo 
a pitchfork bifurcation, reflection symmetry is broken, and traveling waves appear 
{TW curves in Figs. Kllb -d). These traveling- wave solutions manifest themselves as 
fixed points of the three-dimensional ODE system involving Ri, R2 and (/>, but are 
seen to be traveling waves by the fact that the individual phase velocities are nonzero: 
4>i/4'2 = n/m. Since the phase velocity of these waves goes to at the bifurcation, 
they are often called drift waves. We do not consider the stability properties of the 
traveling-wave solutions in this paper. 

The stability results for S'i_2 and S± described above are illustrated in Figs. 13. 213. 41 
for m : n = 2 : 5 and the indicated parameter values. They all satisfy 

s = s' = l, ss'-pp'>0, w' <0, (3.12) 

so that the pure modes S'i_2 bifurcate supercritically in all cases. Next to these plots 
we sketch the type of bifurcation diagram one obtains when increasing p at con- 
stant k along the thin dashed vertical lines. Since in all cases ss' — pp' > 0, both 
mixed modes S± arc stable to amplitude perturbations (sec the explanation follow- 
ing Eqs. H3.8l3.9|l 'l. The stability of S± with regard to phase perturbations depends, 
however, on the eigenvalue e± given by H3.11(l . Because uiy' < 0, this eigenvalue may 
change sign, causing the mixed modes S± to undergo a symmetry-breaking bifurca- 
tion to traveling waves. Figs. 13.213.41 present six different cases characterized by the 
following quantities: 

2 

X^l'-4Ap(5-{-fS') and A = ^ - A^. (3.13) 
\ m / 4o 

The parameter x controls the intersection of the parabolas a = and (3 = 0; they 
intersect if x > and not otherwise. The quantity A determines the relative position 
(//-value) of the minima of the curves a = and f3 = 0. It thus indicates which of 
the two modes, 5*1 or S2, is excited first. S2 appears first when A > 0, while takes 
priority for A < 0; if A = both pure modes onset simultaneously. 

The degenerate case x = A = 0, i.e 7 = Ap = 0, is illustrated in Fig. 13.21 In 
this case the curves a = and (3 = intersect only once and their minima coincide. 
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Fig. 3.2. Stability regions of the pure modes and mixed modes for the resonance 2 : 5 for 
Afi = and 7 = 0. Figures on the right are sketches of bifurcation diagrams associated with the 
vertical paths (dashed lines) in the stability regions on the left. Solid lines correspond to stable states 
and dashed lines to unstable states. T stands for trivial state, (a) 5 = S' = I, s = s' = I, p = 0.4, 
p' = 0.67 u = 0.62 and u' = -1.02. (b) 5 = 5, 5' = 0.5, s = s' = 1, p = 1.5, p' = 0.5 and 
V = —v' = 0.05. In (a) the traveling-wave branch that arises at TW~^ is not shown since we do not 
consider stability properties of traveling solutions. 

While in FiglOi the neutral curve for mode Ai is wider than that for A2, it is the 
other way around in Fig l3.2b . Note that the wavenumber nk/m of A2 is larger than 
that of Ai ; therefore to make the neutral curve of A2 wider than that of Ai requires a 
large ratio of S/S' . Depending on the nonlinear coefficients, all four branches of pure 
modes and mixed modes {Si, S2, S±), or just the pure modes (5'i,5'2) may arise at 
the intersection point of the neutral curves. Because we are considering s = s' ~ 1, Si 
and 52 bifurcate supercritically. For fc = the eigenvalues X^^'^^ (cf. Eqs. 1)3.613.7(1 ) 
determining the stability of Si and S'2 take the simpler form: 

Li-.X^l^ ^fi{l-p')±{iy'\tir^H for5i, (3.14) 

L2 : A^'^ = fi{l - p) for ^2. 
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Near onset (i.e., < /i ^ 1) 5i and 5*2 arc stable if 1 — p' < and 1 — p < 0, 
respectively. In the case m = 2, the stability of Si can be modified at larger values 
of fi by the resonant terms. Near onset the amplitudes of the mixed modes are 
approximated by 

S±: i?? = (m(1 - p) T (i^ - J^'pMm')) (3.15) 
1 — pp' 

Rl = {p{l p') , 

1 — pp' 

Note that small-amplitude mixed modes S± can therefore exist only if (\ — p)(l — p') > 
0. As discussed above, when pp' < 1, the stability of S± is controlled by the phase 
eigenvalue e± of Eq. H3.11|l . Upon substituting Eqs. (|3.15() into Ec^. (|3.11l) wc find 
that 

sign(e±) = sign (t— ^ {niy{l - p') + mv'{l - p))\ . (3.16) 

The results for A/i = 7 = and s = s' = 1 may now be summarized: 

1. If 1 — p > and 1 — p' > both pure and mixed modes bifurcate from 
the trivial state as the forcing p, is increased. The two pure modes are both 
unstable, while one of the mixed modes, either S+ or 5'_, is stable. This case 
is illustrated in the diagram of Fig. 13.2b : in this example sign(e±) = Tl (see 
Eq. H3.16|l ). implying that S+ is stable and S- unstable at onset. 

2. If 1 — p < or 1 — p' < only the two pure modes are present at onset. When 
1 — pp' > 5*1 is stable if 1 — p' < and 5*2 is stable if 1 — p < 0. In particular, 
a bistable situation is permitted. In the example shown in Fig. l3.2l onlv S2 is 
stable because 1 — p' > 0. 

Fig. 13.31 presents two possible unfoldings of the degenerate diagrams shown in 
Fig. 13.21 In both cases the curves a = and /? = intersect twice because x > 0- In 
Fig. 13. 3h , the pure mode 5*2 appears first (A > 0) and is stable. With increasing p, 
however, perturbations in the direction of the other mode, Ai, become increasingly 
important, eventually destabilizing S2 at L2 in favor of the mixed modes S±. Since 
Rl <^ i?2 in the vicinity of L2 and v > 0^ S+ is the stable mixed mode (cf. Eq. H3.11|l ). 
In Fig. 13.3b we have A < 0, and the roles of and S2 are switched. The mode 5"!, 
which is now stable at onset, is destabilized by a perturbation in the direction of A2 
at Lj", generating the mixed state It undergoes a second bifurcation involving 
a perturbation in the direction of A2 at L'^ , leading to the mixed mode Sj^. In this 
case S- is stable, not S+. In contrast with Fig. 13.3b .. however, the pure mode S2 is 
ultimately stabilized at large p since p > s' = 1 . Therefore the cross-interaction term 
p|^2p = p(p + ^m)/s' in Eas. (|2.4l2.5|l . dominates the linear growth rate p of Ai for 
large p and suppresses the perturbations in the direction of Ai. 

Fig. 13.41 shows diagrams for the case where the curves a = and = do 
not cross [x < 0). In Fig. 13.4b the pure mode S2 bifurcates first, while Si does so 
in Fig. 13.4b . The subsequent bifurcations that these states undergo are basically of 
the same type as those in Fig. 13.31 For example, in Fig. 13.4b the mode iS'2 becomes 
unstable to the mixed modes at L2, while Si eventually gains stability at large p 
because p > s'. A feature of Fig. 13.4b that is not present in the previous diagrams is 
the saddle-node bifurcation on the S-i- branch. The appearance of this bifurcation is 
not specific to the case x < but depends on the nonlinear coefficients, in particular 
on the resonance coefficients u and i^' (sec Eq. (|3.5|l '). 
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Fig. 3.3. Stability regions of the pure modes and mixed modes for the resonance 2 : 5. Figures 
on the right are bifurcation diagrams associated with the vertical paths (dashed lines) in the stability 
diagrams on the left. T stands for trivial state Ai = A2 = 0. (a) Same parameters as in Fig. \3.'A i 
but Afi = 0.366 and 7 = 0.7. (b) Same parameters as in Fia. \3./i} but Afj, = —0.5 and 7 = 0.5. The 
branch of traveling waves that arises from TW~ in (a), or that (possibly) connects the S+ with the 
S— solutions in (b) is not shown. 



A comment regarding the validity of the phase diagrams is in order. For the weak 
resonances we are considering in this paper the resonant terms proportional to v and 
v' are of higher order in the amplitudes. Consequently, they have to be considered 
as perturbation terms. Since they are responsible for the splitting of the lines Li and 
TW into and TW^, only those aspects of the results presented in Figures [3. 213. 41 
that persist as this splitting becomes small are expected to hold systematically. Thus, 
the transition to traveling waves at fc = in Fig l3.3b is robust, while that at fc = 
in Fig l3.3h is shifted to ever larger values of /i as the resonant terms become weaker. 
Similarly, the saddle-node bifurcation of 5+ in Fig l3.4b is not robust. It disappears 
through a sequence of bifurcations that involves a merging of the branch 5*+ with Si . 
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For very small |;^| and the branch 5*+ merges with 5*1 close to L 
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Fig. 3.4. Stability regions of the pure modes and mixed modes for the resonance 2 : 5. Figures 
on the right are sketches of bifurcation diagrams associated with the vertical paths (dashed lines) in 
the stability diagram on the left. T stands for trivial state A\ = A2 = 0. (a) Afi = 0.5, 7 = 0.5, 
5 = 5, 5' = 0.5, s = s' = I, p = 0.5, p' = 1.5 and v = -v' = 0.085 (b) Ap = -0.1125, 7 = 0.5, 
S = S' = I, s = s' = 1, p = 1.5, p' = 0.5 and v = 0.5 and v' = —1.01. The branch of traveling waves 
that (possibly) connects the S+ with the S~ solutions in (a), or that arises from TW~ in (b) is not 
shown. 
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4. Side-Band Instabilities. Wc now turn to the analytical core of this paper, 
the stability of the spatially periodic solutions with respect to side-band instabilities. 
Such instabilities are expected to destroy the periodicity of the solutions and provide a 
possible connection with quasi-periodic patterns. We therefore consider perturbations 
of the periodic solutions in the form 

Ai = [l + ai{X,T))Rie'^^^+^^\ A2 = (1 + a2(^, r))i?2e*""''/'"^''^+"^'\ (4.1) 

where aj{X,T) = {a+{T)e'Q^ + aj(T)e^'Q^) with Q ^ and j = 1,2. Note that 
the perturbation wavenumber Q is measured relative to the deviation wavenumbers 
k and nk/m. The linear stability of the pure modes 5i.2 and the mixed modes S± is 
calculated by inserting the ansatz (|4.1|l into Eqs. H2.4l2.5f) and linearizing in . The 
details of these calculations are given in Appendix A. 

4.1. Pure modes Si 2- The linearized system for the perturbations associated 
with the pure modes 6*1^2 separates into two uncoupled 2 x 2-blocks, which allows the 
stability of each pure mode to be calculated analytically. These two blocks can be 
associated with longwave and shortwave instabilities, respectively. The block corre- 
sponding to the longwave instability contains the eigenvalue related to spatial transla- 
tions (i.e. phase modulations). The other (shortwave) block describes the evolution of 
amplitude perturbations in a direction tranverse to the relevant pure-mode subspace. 
We find that in addition to the instabilities discussed in Section 3 the destabilization 
of the pure modes can occur by longwave (Eckhaus) or shortwave instabilities. 

In the case of 5*1 , a straightforward calculation yields 

El : fi - 5k'^ - -fk - {-f + 2k5f/25 = 0, (4.2) 

Ti : {sl3~ p'a)/s + 5'{nk/mf + {v'^\a/sY'/4,5'{nk/mf]^^^^^ = 0, (4.3) 

where Ei is the Eckhaus curve and Fi is the stability limit associated with shortwave 
perturbations. If the curve Fi is crossed, undergoes a steady-state bifurcation with 
a perturbation wavenumber Q given by 

When m > 2 the bracketed terms in Eqs. (|4.3() and H4.4|l are absent. The destabilizing 
eigenfunction then takes the simple form (ai,a2) oc (0, e"'*^™"^) (sec Eqs. H4.1|l and 
(|A.4|I of Appendix A) and allows a correspondingly simple physical interpretation: 
the wavenumber of the destabilizing mode is the wavenumber at the band center of 
the other mode S2 (see Eqs. (|4.1|l '). This is no longer the case for m = 2 because the 
resonance terms (which are linear in A2) affect the stability of 5*1. One consequence 
of this is that the destabilizing mode is composed of two wavenumbers: fc^ ± Q 
with Q given by Eq. H4.4|l . Moreover, it is possible to have one or more non-zero 
wavenumbers fc*, say, for which the perturbation wavenumber Q vanishes. At these 
points (fc,/x) = {k* ^ fx*) the curve Fi merges with the curves (or L^) describing 
stability under homogeneous perturbations. 
In the case of 5*2 we find 

E2: fi + A/i - ?,5'{nk/mf = 0, (4.5) 
r2 : {s'a - pP) /s' + (7 + 25kf/A5 = 0, (4.6) 
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where E2 is the Eckhaus curve and r2 is the stabihty hmit for shortwave instabih- 
tics. As r2 is crossed the pure mode ^2 undergoes a steady-state bifurcation with 
perturbation wavcnumbcr Q given by 

Q^±\k + j/25\. (4.7) 

The associated eigenfunction is of the form (01,02) oc (e-*('=+*)^,0) (see Eqs. (jUTJ 
and (|A.9p ) and, as in the case of 6*1, points to a simple interpretation: the mode 
that destabihzes S2 hes at the band center of 6*1, i.e., its wavenumber is 7/(2(5) = 7 
(cf. Eqs. (|4.1|) . and Appendix A; see also Fig. 12. 1(1 . This result holds for all resonances 
except m : n = 1 : 2, in which case the linear stability of the pure mode S2 is affected 
by the resonance terms (as Si was when m = 2). 

Stability results for the pure modes 5i.2 are shown in Figs. 14.11 and 14.21 In 
Fig. I4.1b .d and Fig. 14.2b the parameters are as in Fig. I3.3b .b and Fig. 13.4b . re- 
spectively, while in the diagrams of Fig. I4.1b .c and Fig. 14.2b we depart from those 
parameters only in setting 7 = 0. The inclusion of these last three cases allows us to 
gauge how much the side-band instabilities are affected by the detuning from perfect 
resonance. 

Note that due to the large ratio S/d' the detuning of 7 = 0.5 has only a small effect 
on the neutral curve of Ai {a ~ 0). Observe that in Fig. I4.1b .b the pure mode 5*2 
is stable within the region bounded by the curves E2 and r2 while Si is everywhere 
unstable due to shortwave instabilities. In contrast. Fig. I4.1b .d depicts a situation 
with stable regions for both pure modes. Note that the pure mode S2 experiences 
only shortwave instabilities as the forcing is decreased while the pure mode 6*1 can lose 
stability to either longwave or shortwave perturbations, provided the wavenumber k is 
not within the interval defined by the merging points of Fi and (see Fig. 14.1b and 
I4.1t l'): over that interval the instability suffered by S'l is to perturbations preserving 
the periodicity of S'l and gives rise to (steady) mixed modes. The parameters of 
Fig. l4.2l do not allow for codimension-two points (intersection of the curves a = and 
(3 = 0). Despite this difference, the stability results for S2 are qualitatively similar to 
those of Fig. im The effect of 7 on S'l, however, is more striking than in Fig. l4.1b -d. 
In particular. Si can undergo a shortwave steady-state instability (along Fi) when 
7 = 0.5 but not when 7 = 0. 

4.2. Mixed modes S±. The stability of the mixed modes S± is governed by 
four eigenvalues which must be determined numerically. Two of these are associated 
with amplitude modes and are always real. The remaining two may be real or complex 
and are related, respectively, with the translation mode and with the relative phase 
between the two modes Ai and A2. Throughout this stability analysis we focus 
on the influence of the detuning parameter 7 on the (side-band) instabilities of S±. 
Furthermore, because of the invariancc of Eqs. (|2.4I2.5|I under the transformation 

7^-7, X^-X, 

we may take 7 > without any loss of generality. 

The results of this stability analysis are shown in Fig. 14.11 which depicts a situ- 
ation with two codimension-two points, and Fig. 14.21 which presents a case with no 
codimension-two points. Since this difference seems to have only a minor effect on the 
stability of S± we focus our discussion on the case of Fig. im fi.e. two codimension-two 
points). 
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Fig. 4.1. Stability regions of the pure modes 5i,2 o.nd the mixed modes S± for the resonance 
2 : 5 and the following parameter sets: (a) Afi = 0.366, S = S' = I, s = s' = I, p = 0.4, p' = 0.67, 
u = 0.62, u' = -1.02 and 7 = 0; (b) as in (a) but 7 = 0.7 (cf FiaU^): (c) Ap. = -0.5, <5 = 10, 
5' = 0.5, s = s' = 1, p = 1.5, p' = 0.5, u = -u' = 0.05 and 7 = 0; (d) as in (c) but 7 = 0.5 
(cf. Fia iy.yb ). The symbols in b) and d) indicate the parameters used in the numerical simulations 
described in Sec. 

4.2.1. Mixed mode S+. When ^ = fFig. l4.1b .c). the mixed mode .S*^ becomes 
unstable to short-wavelengtli perturbations on the part of r+ closer to L2, while a 
longwave stability analysis captures the part of r+ closer to TW^ . Both instabilities 
are of steady-state type with the one close to TW^ breaking the reflection symmetry 
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Fig. 4.2. Stability regions of the pure modes Si,2 and the mixed modes S± for the resonance 
2 : 5 and the following parameter sets: (a) Afi = 0.5, 7 = 0, (5 = 5, <5' = 0.5, s = s' = I, p = 0.5, 
p' = 1.5, 1/ = -I/' = 0.085. (b) as in (a) but 7 = 0.5 (cf Fig^^\i). 

of the pattern. This is expected to lead to a drift of the pattern; additional calculations 
for other parameter sets (not shown) suggest that this situation is characteristic of 
S-\- and 7 = 0. 

When 7 7^ the mixed mode S+ displays both steady-state and Hopf bifurca- 
tions. The two relevant eigenvalues are associated with the translation mode and 
the symmetry-breaking mode, respectively, with the latter eigenvalue going to zero at 
TW~^. The change in character of the instability along is illustrated in Fig. 14.31 
Figs. I4.3b .c give the growthrates of the dominant modes as r+ is crossed at specific 
points marked in Fig. 14.3b . At points 1-5 the bifurcation is steady but as one continues 
in a counterclockwise direction (points 4-6) the eigenvalues merge and become complex 
over an interval of Q-values. This interval containing complex conjugate eigenvalues 
continues to grow, leading eventually to an oscillatory instability superceding the 
steady one. A codimension-two point therefore exists near {k, ji) = (—0.155, 0.581) at 
the transition from steady-state to Hopf bifurcation; at this point there are two un- 
stable modes, at different Q, one steady and one oscillatory. As one goes still further 
in the counterclockwise direction, the Q-band over which the eigenvalues arc complex 
reaches (5 = and the oscillatory instability becomes a long-wave instability. This 
interaction was studied in the context of Taylor vortex flow 00] and has been found 
in many others problems [381 r20| . 

Similar transitions occur if ji and 7 are varied (rather than /i and k). We illustrate 
this in Fig. 14.41 bv plotting the growth rate of the most unstable perturbations after 
crossing, at fixed k = —0.1, the upper part of r+ which is close to TW^ fFig. 14.4b ). 
and the lower part of r+ which is close to L2 (Fig. 14.4b ). Along the upper part of 
r+ the mixed mode 5'+ typically sees a shortwave oscillatory instability; for small 
7, S+ can also lose stability to longwave perturbations. The lower part of r_|_ is 
characterized by either steady-state or oscillatory instability of shortwave type (see 
Fig. 14. 4b V There is again a codimension-two point, near (7, fi) ~ (1.035, 0.47), where 
Hopf and steady-state bifurcations coalescence. For 7 less (greater) than this critical 
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Fig. 4.3. Growth rates of the two most unstable perturbations of S+ when r+ is crossed at the 
points indicated in the upper figure: Re{Xi) (solid-thin line) and Re{\2) (dashed-thick); Ai and A2 
are complex where there is a single curve, i.e., where Re{\\) = iJe(A2). We use 7 = 0.4 with the 
remaining parameters as in Figs. \4.1}p. 



value the bifurcation is steady (oscillatory). The route to the oscillatory instability is 
as described above: with increasing 7 the two modes associated with translation and 
reflection symmetry interact more and more, ultimately leading to a Hopf bifurcation. 

The steady-state shortwave instability of S+ (see Fig. 14.4b ) always appears in the 
vicinity of L2. This proximity implies that the amplitude Ai of the mixed mode is 
very small. Consequently, the behavior of 5"+ in this region is similar to that of the 
pure mode 82- The growing perturbations are predominantly in the direction of Ai 
and the associated wavenumber Q is given in a first approximation by Eq. (|4.7|l . i.e., 

Q~±\k + -i/25\. (4.8) 

As demonstrated in Fig. 14.51 Eq. (|4.8|) provides a good approximation over the 
interval 0.1 < 7 < 0.7 when |fc| < 0.1. It is interesting to note that, when 7 is not too 
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Fig. 4.4. Growth rate of the most unstable perturbation of S+ when r+ is crossed vertically 
at k = —0.1 for the indicated values 0/7, (a) on the upper part 0/ r_|_ (close to TW^ ) and (b) on 
the lower part (close to L2 )■ The remaining parameters are as in Fig, \4.ll i-b. Solid (dashed) lines 
correspond to real (complex) eigenvalues. 
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Fig. 4.5. Dependence of the perturbation wavenumber Q of the steady mode destabilizing 
on the detuning parameter 7 for the resonace case m : n = 2 : b. Parameters as in Fig \4.1\ i,b with 
a) k = —0.1 and b) k = 0. Open circles correspond to the perturbation wavenumber Q calculated 
numerically from Eq. IB.<HB.3D . Solid lines correspond to the indicated approximations. 
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small, Eq. H4.8|l applies independently of the resonance m : n. The robustness of this 
result is not unexpected since the resonance terms do not influence at first order the 
stability properties of the mixed mode 5*+ (because \Ai \ is small). However, Eq. (|4.8|l 
ceases to apply when 7 becomes very small because resonance effects start to play a 
significant role in the wavenumber selection of the linear response when crossing the 
lower part of r_|_. In particular, for fc = and 7 <C 1 one can show 

Q ~ (7V)i/4 [r2o(2s7(ss' - pp')f/^] /2S + 0(7"+'). (4.9) 

where r2o denotes the amplitude of S2 (with fc = 0) on the curve ^2- The derivation of 
Eq. 1)4. 9() . given in Appendix B, relies on the fact that for 7 = the curve r+ becomes 
tangent to the curve L2 at fc = 0. Fig. 14.61 shows that 1)4. 9|l is in excellent agreement 
with the numerical results obtained directly from the characteristic equation l|B.2|l 
and applies to other weak resonances with m : n ^ 2 : 5. Note that for larger values 
of n Eq. H4.9|l is valid over a wider range of 7 values (since the error is 0(7"+^)) and 
the cross-over to the behavior H4.8|l is shifted to larger values of 7. 




Fig. 4.6. Dependence of Q on the detuning parameter 7 for the indicated resonances m : n 
and fc = 0. Remaining parameters are = 0.366, & = &' = 1, s = s' = 1, p = 0.4, p' = 0.67, 
y = 0.62 u' = —1.02. Open circles correspond to numerical results of Eq. IB. SIB. 31) and lines to the 
approximations \4.<^ and {4.yt , Q oc 7/25 (dashed) and Q oc 7"/'* (solid), respectively. 



4.2.2. Mixed mode S_. As illustrated in Fig. 14.1b . the mixed mode S- can 
be everywhere unstable when 7 = 0. For the parameters of Fig. 14. lb . such a situation 
persists until 7 w 0.6, when a stability region for 5_ emerges. This stability region 
widens as 7 is increased (see Fig. 14.1b '). 

As in the case of S'+, the two most important eigenvalues for the stability of S- are 
those associated with the translation and parity-breaking modes. The oscillatory in- 
stability results from an interaction between these two eigenvalues and, consequently, 
is found in the vicinity of the bifurcation set TW~ ; the steady (shortwave) instability 
occurs near the initial bifurcation producing the mixed mode 5_ from the pure mode 
(i.e., near L^). Thus, the mixed mode 5_ undergoes the same kinds of transitions 
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described above for 5*+. Depending on where r_ is crossed in the (k, /i)- or the (7, /i)- 
plane, the instabihty may be cither osciUatory or steady, and either of fong-wave or 
short-wave type. 




Fig. 4.7. Growth rate of the most unstable perturbation when F- is crossed vertically at k = 0.5 
for the indicated values of ^ and the remaining parametersas in Fig, \4.ll i-b. (a) on the upper part 
of T— (close to TW~ ) and (b) on the lower part (close to ). A solid line is used when the 
associated eigenvalue is real and a dashed line when it is complex. 



Note that on the upper part of r_ in Fig l4.7b . as 7 is increased the two eigen- 
values merge and become complex at ever smaller values of Q. In contrast, on the 
lower part of r_ raising 7 increases the damping of the parity-breaking mode, thus 
discouraging its interaction with the translation mode, this behavior is opposite to 
what occurs with 5+. It is also interesting to note that variations in the detuning 
parameter 7 do not induce a significant change in the perturbation wavenumber Q 
associated with the bifurcation. This is because the instability of S- takes place close 
to where the destabilizing perturbations are (preferentially) in the direction of 
A2. The perturbation wavenumber Q associated with the instability is given, to first 
approximation, by Eq. (|4.4() . Thus Q depends mainly on the resonance ra : n and 
only implicitly on the parameter 7 (through the function a defined in Eq. Ij^.^jl 'l. 

5. Numerical Simulations. The main objective in this section is to investigate 
the nonlinear evolution of the side-band instabilities of the pure and mixed modes 
discussed in Section 01 The results are obtained for a periodic domain of length L 
using a finite-difference code with Crank-Nicholson time-stepping. Each simulation 
of H2.4I2.5|) uses as an initial condition an unstable periodic steady state (i.e., a pure 
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or mixed mode of the type discussed in Section |2Jl with a small random pertm'bation. 
The final solutions obtained in this way are typically 'quasi-periodic' although periodic 
solutions are also found. Here we use the term quasi-periodic to refer to any solution 
for which the wavenunibers of Ai and A2 are not in the ratio m/n. Therefore the full 
reconstructed field H2.2(l will not be composed of wavenumbers in the ratio m/n; in 
general it will be quasi-periodic. 

First we describe the evolution of system H2.4I2.5() when the pure modes become 
unstable. If the instability is longwave it affects only the excited amplitude and the 
solution remains in the pure-mode subspace. For example, in the case of Si (see 
Fig. [STj the Ion gwave instability involves Ai but not A2, which remains zero. As 
expected, the changes affect predominantly the phase of Ai and evolve into a series of 
phase-slips. The final state, which is time-independent, is shown on the right part of 
Fig. 15.11 In the case of 5*2 (not shown) the longwave instability behaves analogously: 
it is now A2 that changes, while Ai remains zero. 




Fig. 5.1. Evolution of the longwave instability of the pure mode Si. Space-time diagram of the 
real part of Ai (left) and final state (right) with thick (thin) lines denoting the real (imaginary) part 
of Ai. Parameters as in Fig. \4.l[ i with k\ = k = 0.15, fj, = 0.37 (indicated by a square in Fig, \4.l\ l) 
and L = 250. 



Depending on the detuning parameter 7 the transition from the pure modes to the 
mixed modes can also involve side-band instabilities. For the parameters correspond- 
ing to the stability diagram shown in Fig. I4.2l the pure mode Si becomes unstable to 
the mixed mode S- via a side-band instability for 7 7^ and k > 0.15 (cf. Fig. 14.2b '). 
while for 7 = 0, i.e. when the two critical wavenumbers are resonant, the instability 
preserves the periodicity of the solution for all k (cf. Fig. 14.2b '). The final states 
that are reached after a small random perturbation has been applied to a pure-mode 
solution 5*1 are depicted in Fig l5.2l The movies fig5. 2a.avi and fig5. 2b.avi show 
the corresponding temporal evolution of A2 . The sequence of phase slips that is ev- 
idenced in the sharp peaks of the local wavenumber during the early stages of the 
evolution is a consequence of the random initial perturbation. After this transient the 
relevant perturbation mode starts to dominate. In the case 7 = 0.5 that mode has a 
strongly modulated wavenumber reflecting the fact that the wavenumbers of Ai and 
A2 do not satisfy the resonance condition m : n = 2 : 5, resulting in a quasi-periodic 
reconstructed solution u{x) (see also discussion leading to (15. 1|) below). 
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Fig. 5.2. Final state after a small random perturbation of the pure mode Si. Parameters in 
(a) and (b) as in Fig. \4.i^ ,b, respectively with L = 250 and 7 = 0, A; = 0.1, fi = 1.4 in a) and 
7 = 0.5, k = 0.2, fi = 1.4 in b). The corresponding movies [movie:fig5. 2a.avi] and [movie:fi,g5. 2b.avi] 
shovj the temporal evolution of A2. Red and yellow lines show the real and imaginary part of A2, 
respectively, white its magnitude \A\, and green the local wavenumber. 



The longwave instability of the mixed modes affects Ai as well as A2 (see Fig. l5.3f) . 
As with the pure modes, the phase perturbations evolve into phase slips, which change 
the wavenumber and lead eventually to a stationary state. In the case shown in 
Fig. 15. bi phase slips occur only in mode ^2- As in the mixed mode shown in Fig l5.2b . 
the strong deviations of ^1^2(2;) from purely sinusoidal behavior, which are due to the 
resonant terms proportional to ly and i/', indicate that the reconstructed solution u{x) 
is not periodic. 

The mixed modes can undergo a shortwave steady instability as well. As discussed 
in Section^ it occurs when the stability limit F^ (F_) of S+ {S-) is close to the 
transition line L2 (^T)- T^^^ instability is primarily in the direction of the mode 
with smaller amplitude: Ai for r_|_ near L2 and A2 for F_ near L'^ . The resulting 
evolution of system H2.4I2.5|) is shown in Fig. l5.4l for parameters near F^. It confirms 
the expectation that at least initially only Ai changes, but not ^2- The space-time 
diagram for Ke{Ai{X, T)) in Fig. l5.4l shows how the growing perturbation determines 
the wavenumber of the final state; this change from the initial wavenumber takes place 
via phase-slips. 

The wavenumber fci of Ai in the final state is controlled by the detuning parameter 
7 in the manner predicted by the linear analysis: fci ~ j/(26). In other words, the 
number of maxima observed in Ai is A^i ~ Lj/{Att5) (see lower figures of Fig. I^T^ . 

The evolution of Ai and A2 after a shortwave steady instability on r_ near Lj~ 
is shown in Fig. 15.51 This time the amplitude Ai remains nearly constant while 
A2 changes. The destabilizing mode was shown in Section ^ to be of the form 
aJ(T)e*('="+'3)-^ + a2{T)e'^''^-'^'>^ with Q defined by Eq. g^. For the cases il- 
lustrated in Fig. 15.51 we have |a^(T)| <^ |a^(r)| and the wavenumber of the linearly 
unstable mode is approximately given by k2 ~ fc-^ ~ IQI- Note that this wavenumber 
is essentially determined by the resonance terms and would be without them (see 
Eg. I4.4|l . In particular, for the example considered in Fig. 15.5b . where ~ 0.62 and 
u' = -1.02, we have |fc2| ^ 0.144, while in Fig.lCT) u = -ly' ^ 0.05 and |fc2| ~ 0.018; 
note that with regard to this effect both cases are are analogous, except for the mag- 
nitude of the resonance coefficients. 

In addition to the wavenumber ^2 determined by the linear stability analysis. 
Fig. 15.51 reveals a second prominent wavenumber fcg, which is the result of the res- 
onance term A™~^A^. It acts as a driving term for A2 and generates a mode with 

23 



SB 




0.3 0.4 0.0 o.a I 

^ Fina] statp 

i ^ 

a 



02 



0.4 OB' 



S 1 

I 





— 1 — 


1 


1 




1 






>^^^ 


D D.I 


D.Z 


D.3 


D4 


D5 


D.6 


D.7 


DB 


D.3 1 





0.1 



D.Z D3 D.4 D5 



D6 D.T 



D.B 0.3 1 
JE/i 



Fig. 5.3. Evolution of longwave instability of the mixed mode S+. Shown are space-time 
diagrams of the real part of the amplitudes and the final state with thick (thin) lines denoting the 
real (imaginary) parts of Ai and A2, respectively. Parameters as in Fig.\J^ with k = 0.1, fi = 0.423 
(indicated by rhomb) and L = 300. 



wavenumbcr that is determined by 

nki ± (m — 1)1^2! 



±\k'\ 



(5.1) 



where ki is the wavenumber associated with amphtude Ai. In the evolution shown 
in Fig l5.5l ki remains unchanged. The same mechanism is at work for the short- 
wave instabihty of S+ (cf. Fig l5.4|l . which occurs near L2. There the wavenum- 
ber modulation of Ai is, however, negligible. It is driven by the resonance term 
A'^~^A'!p, which is proportional to the fourth power of the amplitude Ai [n = 5) 
which is very small where S+ branches off the pure mode ^2. This is not the case 
for the short-wave instability of 5_ fFig l5.5() . which occurs near Lj~. There it is 



A2 that is small; but it enters the resonance term A™" A" linearly (m = 2) and 
therefore the resonance term provides a relatively strong modulation of A2. Sub- 
stituting the final values of ki (same as the initial value) and fc2 (calculated above) 
into Ea. (|5.1|) one obtains jfcjl ~ 2.644 in the case of Fig. 15.5b and |fc2| = 0.63 for 
Fig l5.5b . For the numbers (A^i,iV2,A^2) of wavelengths associated with fci, fc2, and 
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k'2, respectively this implies (TVi, iVj, 7V^) = (20,6,106) ~ L/(27r)(0.5, 0.144, 2.644) for 
Fig. ESli (where L = 250) and {Ni,N2,N^) = (6,1,31) ~ L/(27r)(0.124, 0.02, 0.63) 
for Fig. 15., (where L = 305). These results eompare quite well with the results in 
Fig.ESlwhere (TVi, iVz, 7V^) = (20,6,106) in Fig.ESt and (iVi, iVz, 7V^) = (6,0,30) in 
Fig.EHb. 
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Fig. 5.5. Space-time diagram showing the real parts of Ai and A2 when mixed mode S- 
undergoes a shortwave steady instability. The parameters in (a) are those of Fig. \4-l\ > (circle) with 
k = 0.5, /I = 1.05 and L = 250. Parameters in (b) are as in Fig. (circlejwith k = 0.125, 

^l = 0.95 and L = 305. 
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Qualitatively different is the evolution of the oscillatory instability of the mixed 
mode, as illustrated in Figs. 15. 615. 91 Small perturbations initially evolve into standing- 
wave oscillations. They do not last for long, however, and decompose into left- and 
right-traveling disturbances, which then form localized low-wavenumbcr domains of 
waves drifting to the right and to the left, respectively. Such localized, propagating re- 
gions of 'drift waves' are typical for parity-breaking instabilities because the extended 
drift waves are generically unstable at onset |4()l El ^] and may become stable only 
for larger amplitudes 0]. The localized waves can be described using equations for 
the amplitude of the parity-breaking mode and the phase of the underlying pattern 
and can arise stably when the parity-breaking bifurcation is subcritical j^l 12 1| as well 
as when it is supercritical |4U[ IS) . They are related to the solitary modes observed 
in experiments in directional solidification [431 117| as well as in Taylor vortex flow 
|4b| . and in simulations of premixed flames 0]. In these systems the parity-breaking 
bifurcation arises from a 1:2 mode interaction |26l I4(J| rather than the 2:5-resonance 
considered here. 

The localized drift waves do not always persist. In the simulation shown in Fig l5.6l 
a sequence of phase slips occurs in mode A2 at the trailing end of the localized drift 
wave. The phase slips substantially reduces the wavenumber in the growing domain 
between the location where the localized drift wave first is created (at x « 0.6/L) 
and its trailing edge. Since no phase slips occur in mode Ai the wavenumbers of 
Ai and A2 are not resonant any more in this growing domain, which results in a 
strong modulation of the wavenumber. Due to the periodic boundary conditions the 
growing, localized drift wave collides eventually with this domain and is absorbed by 
it. After that the pattern becomes stationary. The stationary domains with strongly 
modulated wavenumber do not always absorb the localized drift waves. Fig l5.7l and 
the movie f ig5 . 7 . avi show a case in which the localized drift waves at times are also 
reflected by the stationary domains or pass through them. We have not investigated 
which factors determine the outcome of the collisions. In the simulation shown in 
Fig l5.7l the system eventually evolves into a stationary state. It consists of relatively 
large domains in which the wavenumber ratio is very close to 2:5 and small domains 
in which the wavenumbers are strongly modulated and not in this rational ratio. 
One may have expected that the localized stationary domains attract each other and 
eventually merge; but this was not observed. Presumably, the strong wavenumber 
oscillations, which are particularly visible in A2, lock the domains in place. Thus, 
in the final, stationary state shown in Fig l5.7b one has a periodic pattern that is 
interrupted by small domains in which the pattern is not periodic. 

The transients and the resulting final state depend very sensitively on the am- 
plitude of the initial perturbations of the periodic state. The relatively large pertur- 
bations used in Fig l5.7l lead to a large number of localized drift waves, which upon 
their collision generate localized stationary structures. Since the localized stationary 
structures often absorb the localized drift waves the final state tends to be stationary 
in this case. If the initial perturbations are much smaller, fewer drift waves arise. 
If they travel in the same direction no stationary localized structures are generated 
and the localized drift waves persist. Such a case is illustrated in Fig l5.8l5"^ and 
movie fig5.8.avi. Here, uniformly distributed perturbations with maximal ampli- 
tude 0.001 were applied independently to the real and imaginary part of both ampli- 
tudes. In this case only a single domain forms. It develops a defect in its interior, 
which eventually disappears at its trailing end (cf. movie fig5.8.avi). Initially, 
the propagation velocities of the leading and the trailing front of the domain are not 
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equal and the domain grows. As time goes on, however, the front velocities converge 
to a common value and one sees a stable, localized, propagating domain of fixed size 
containing traveling waves (see Fig. I5.9|l . It is interesting to note that the wavenum- 
bers selected in these traveling domains are in 2 : 5 resonance, while the stationary 
regions have slightly shifted wavenumbers that are not in resonance (see upper panel 
of Fig. 1^7^ . The physical solution (cl. ^2.'2^ ) can therefore by described as consisting of 
regions of localized periodic waves traveling through a stationary quasi-periodic pat- 
tern. This type of solution is generally found when the side-band oscillatory instability 
occurs very close to the parity-breaking bifurcation. To describe the parity-breaking 
bifurcation in more detail the appropriate amplitude-phase equations [^1 1401 121 0] 
could be derived from the coupled Ginzburg-Landau equations. They would, in par- 
ticular, allow a better understading of the localized drift waves. Such an undertaking 
is, however, beyond the scope of this paper. 

Finally, in order to determine the character of the bifurcations we have also inte- 
grated Eqs. H2.4I I2.5|l using the final states as initial conditions and varying /i so as 
to take the system back into the stable regions of mixed modes. The results obtained 
indicate that the instabilities discussed in this section are subcritical (except for the 
cases shown in Fig l5.2|l . The solutions persist over a large interval of the forcing 
giving rise to bi-stability between periodic and quasi-periodic states. This behaviour 
has been found in all cases we have considered. Therefore we expect that the coexis- 
tence of periodic and quasi-periodic stable states may be regarded as a characteristic 
property of the system Eqs. ()2.4I I2.5|l . 
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Fig. 5.6. Space-time diagram of the phase- gradient and the magnitude of the amplitudes Ai 
and A2. Lower figures: Real (thick) and imaginary (thin) part of Ai and A2 corresponding to the 
final stable state. Parameter as in Fig. \4-ll > (triangle) with k = 0.5, /i = 1.308 and L = 250. In the 
gray-scale plots dark (light) stands for low (high) amplitudes. 
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Fig. 5.7. aj Space-time diagram of the phase- gradients amplitudes Ai and A2, respectively. 
Parameters as in Fig. \4.1p but 7 = 1 and with k = 0.5, fi = 1.42 and L = 250. Lower figures 
correspond to the final state: a) Real and imaginary part of the amplitudes Ai and A2 b) phase 
gradients of Ai and A2 . Dotted line indicates perfect 2 : 5 resonance. The temporal evolution of A2 
is also shown in movie [movie: fig. 5. 7.avi]; red and yellow lines show the real and imaginary part of 
A2, respectively, white its magnitude \A\, and green the local wavenumber. 
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Fig. 5.8. Space-time diagram for the magnitude and the phase-gradient of the amplitudes Ai 
and A'2,. Parameters: as in Fig. |4- -^^ ^'^t 7 = 1.3 and with k = 0.5, = 1.42 and L = 250. In the 
gray-scale plots dark (light) stands for low (high) amplitudes. The temporal evolution of A2 is also 
shown in movie [movie: fig. 5. 8.avi]; red and yellow lines show the real and imaginary part of A2, 
respectively, white its magnitude \A\, and green the local wavenumber. 
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6. Conclusions. In this paper we have investigated a near-resonant steady m : n 
mode interaction in one-dimensional dissipative systems with reflection and transla- 
tion symmetry using coupled Ginzburg-Landau equations that describe the slow evo- 
lution of the two relevant mode envelopes. One of the goals was to shed some light on 
the competition between periodic and quasi-periodic patterns by addressing in detail 
the stability of the periodic solutions, emphasizing particular instabilities that can 
lead to quasi-periodic patterns. 

The simplest stationary solutions of this system are reflection-symmetric periodic 
patterns: the pure modes Si^2 and the mixed modes S±. Building on the results 
of Dangelmayr ^3] we established the stability properties of these solutions with 
respect to perturbations that preserve the periodicity of the solution. Among the 
many possible choices of parameters we restricted our attention to sets of values that 
produced two supercritical primary bifurcations to pure modes followed by secondary 
bifurcations to mixed modes - these mixed modes, in turn, suffered symmetry-breaking 
bifurcations to traveling waves. 

Side-band instabilities modify the stability regions of the pure and mixed modes. 
In particular, we find that these steady periodic patterns can undergo longwave as well 
as shortwave instabilities. For the pure modes the side-band instabilities are always of 
steady-state type. The destabilizing longwave perturbations correspond to the well- 
known Eckhaus instability and depend only on the linear coefficients of Eqs. (|2.4I2.5|I . 
In contrast, the shortwave instabilities depend on nonlinear interactions and on the 
resonance considered (i.e., on m and n). In particular, we find: 

1. If 771 > 2 and a pure mode, 5*1 say, loses stability in a shortwave steady-state 
bifurcation, the wavenumber of the growing perturbation lies at the band center of 
the other mode (i.e. of 52)); this instability is affected by the detuning parameter 
7 = 7/2(5, which measures the mismatch between the two critical wavenumbers. 

2. For wavenumber ratios of the form 2 : n one of the resonant terms enters 
into the stability calculation for the pure mode ^i. Away from the band center this 
pure mode can experience shortwave instabilities while close to the band center it 
can lose stability only to homogeneous or longwave perturbations. In the short-wave 
case, due to the resonant terms the wavenumber selected by the instability is shifted 
with respect to the band center of the other pure mode, ^2. As the band center is 
approached the perturbation wavenumber of the short-wave instability can go to 0. 

Close to their bifurcation from the pure modes, the mixed modes inherit the 
stability properties of the pure modes with respect to short-wave perturbations. In 
addition, they may undergo an oscillatory instability. It results from the interaction 
of the translation and the symmetry-breaking modes and can be strongly affected by 
the wavenumber mismatch. 

The nonlinear evolution ensuing from the instabilities was investigated through 
direct numerical simulation of Eqs. H2. 412. 5(1 with periodic boundary conditions for the 
slowly varying envelopes. It was found that the transition from the pure modes to the 
mixed mode can depend strongly on the wavenumber mismatch. We have identified 
cases in which without a wavenumber mismatch the instability of the pure mode leads 
to a periodic mixed mode, while in the non-resonant case the instability of the pure 
mode is determined by a side-band instability resulting in a quasi-periodic pattern. 

Side-band instabilities of the mixed modes generically destroy the spatial periodic- 
ity of the stationary, reflection-symmetric states and produce quasi-periodic patterns. 
They are typically characterized by a spatial modulation of the local wavenumber with 
a wavelength that is determined by the detuning parameter (in a manner predicted 
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by the linear analysis) and the resonant terms (as indicated in Ea. (|5.1|l 'l. The size 
of the modulation depends strongly on the strength of the resonance terms. For the 
weakly nonlinear 2:5 resonance discussed here the effect is only noticeable for 5_ (cf. 
Fig l5.5() : for S+ the resonant terms are too small (cf. Fig l5.4|l . 

The oscillatory instability tends to lead to the formation of propagating fronts 
separating reflection symmetric patterns from drift waves with broken reflection sym- 
metry. The fronts may interact to form stable localized drift waves. These may sub- 
sequently collide with each other, destroying the initial periodic pattern completely 
in some cases and only partially in others. In the latter case the system relaxes to a 
very interesting pattern composed of alternating localized stationary domains of pe- 
riodic and quasi-periodic states. We find other cases where all the fronts travel in the 
same direction, eventually achieving the same velocity and producing stable localized 
propagating domains of drift waves. The wavenumbers in these traveling domains 
are selected by the resonance condition (i.e., they are in the ratio m/n), whereas the 
wavenumbers associated with the surrounding steady periodic pattern are no longer 
resonant. Thus, the resulting patterns can be described as localized periodic traveling 
wave structures propagating through a stationary quasi-periodic background. 

In this paper we have focused on the periodic solutions and their stability. Com- 
plementary to this analysis would be a study of the quasi-periodic solutions. Of 
course, the set of solutions in which the two modes have non-resonant wavenumbers 
is considerably larger than that of the periodic solutions and such an analysis is be- 
yond the scope of the present work. In preliminary numerical simulations we have 
considered the stability of certain quasi-periodic solutions and found that the range in 
the wavenumber mismatch over which they are stable is smaller than that of periodic 
solutions with nearby wavenumbers |22j . Thus, at least in the regime considered, the 
resonant terms enhance the stability of the periodic solutions. Another promising 
line of investigation would be the reduction of the coupled Ginzburg-Landau equa- 
tions H2. 412. 5(1 to two coupled phase equations. This should be possible for long-wave 
perturbations of the mixed-mode solution if the smallness of the resonant terms is 
exploited and it is assumed that such terms are of the same order as the gradients in 
the magnitude of the amplitudes (cf. (|2. 612. 7(1 ). 

We gratefully acknowledge discussions with J. Porter and P. Umbanhowar. 

Appendix A. Stability of Pure and Mixed Modes. 

In this appendix we provide some of the details of the stability analysis of the 
pure and mixed modes, presented in Sec^ 

A.l. Pure modes. The perturbed pure mode Si is written as 

Ai = i?ie'^^ (1 + a+(T)e*^^ + a^iT)e-"^^) (A.l) 
A2 = e'''^^ {a+(T)e'^^ + a^(T)e-''5-'^) . (A.2) 

Inserting this ansatz in ((2.412.5(1 and linearizing in and we obtain 

a± = (a - si?2 - SQ'^ T (7 + 26k)Q) a± - sRlaf, (A.3) 
d± = (/3 - p'Rj - T 2S'Q{nk/m)) + {i^'R^a^}^^^ , (A.4) 

where the bracketed term with subscript m = 2 is present only if m = 2, and a and 
P are defined by Eq. ((3.3(1 . The eigenvalues associated with Eq. ((A. 3(1 are 



A± = -{a + 5Q'-) ± v/a'+Q2(^ + 2fcJ)2, 
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(A.5) 



and with Eq. HA.4p are 

Eqs.ljXSj and lfO|) are related to longwave and shortwave instabilities, respectively. 
Note that for m > 3 the perturbation amplitudes and decouple and the desta- 
bilizing mode is composed of just one wave of the form (a^, 0) or (0, a^). In contrast, 
due to the resonant term when m = 2, the eigevectors have (in principle) components 
in both directions and the destabilizing mode is composed of two waves with different 
wavenumbcrs. 

For the pure mode 5*2 we have 

Ai = e'''^ (a+ (r)e"5^ + {T)e-'^^) (A.7) 
A2 - i?2 e'^^^ (1 + a+{T)e''^^ + a^{T)e-'Q^) , (A.8) 

where af and satisfy 

af = {a- pR^-SQ^ Til + 2Sk)Q)af, (A.9) 

af {l3- 2s' Rl - T 25'{nk/m)Q) - s'Rjaf. (A.IO) 

The stability of 6*2 is determined by two eigenvalues: 

\f = {a-p^~SQ')±\Q\\j + 2k6l 



A± =-(/? + 5'02)±y^2+ 4^/2 (£fc)'Q2^ (A.ll) 

which are associated with the shortwave and longwave instabilities, respectively. 

A. 2. Mixed modes. In this case we consider the system ()2.6|) - (I2.9|) for the 

evolution of the real amplitudes Rj > and phases (for j = 1,2). Thus (abusing 
notation) we write 

RiiX,T)^Ri+ri{X,T), MX,T) ^ kX + ^^{X,T) (A.12) 



where 



R2{X,T) = R2+r2{X,T), MX,T) ^ -kX + ip2{X,T), (A.13) 

m 



r, = [r+iT) - zr-(T)] e^^^ + [r+(r) + zr-(T)] e 
VP, = [^+(T) ~ ^^-iT)] e^Q^ + [^^{T) + ^^JiT)] 



-iQX 



for j = 1,2. The linearized problem for the perturbations and is given by 



/ ai 

-2SQ{k + 7) 

C2 

V 



-2SQ{k + -f) 
bi 

d2 
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where 



ai = - Sk^ - -fk - SQ^ - (3si?? + pRl) + v{n - 1)R'^~^R^ cos(n0i - mi^a), 

as = -2pi?ii?2 + i'mR'l'^R^'~\ 

61 = -vnRl-'^R'i' - 5Q^, 62 = vmRl-^R'^-\ 

ci= n-Afj.- 5'{kn/mf - - (Ss'RI + p'Rj) + v'{m ~ l)R'^-'^R'^ cos(n(^i - m^z), 
C2 = -2p'i?ii?2 + v'nR'l-^R'^^'^, 
di = -v'mR'lR'i'-^ - (5'Q^ ^2 = i^'nR'l~^ R^'-\ 
Appendix B. Derivation of (|4.9|> . 

To derive H4.9|l we first consider the characteristic equation in the form 

+ G-iX^ + G2A2 + GiA + Go = 0, (B.l) 

where Gj (j = 0, . . . , 3) arc polynomials in whose coefficients depend on the 
parameters directly and also through the amplitudes of the mixed modes. When the 
bifurcation is steady and of shortwave type the location of r_|_ and the wavcnumbcr 
Q ^ are obtained as solutions of the system Go = dGo/dQ^ = 0, which yields 

(5152)' - 4.90.92 - _ 18go5i52 = (B.2) 

for the location of the lower part of r+ and 

= (5152 - 95o)/(65i - 252) (B.3) 

for the perturbation wavenumber of the most unstable perturbation. Here the funtions 
go, gi and g2 can be written as 

go = -{aici - a2C2){S'bi - Sdi) - 7^cidi, 

gi = 5S'{aiCi — 0202 + airfi + ci6i) + 5'^cidi + 5'^aibi + 7^(S'(ci + di), 
32 = -S^S'ici + di) - 6'^'d{ai + bi) - -f^S'\ 

where a^, bj, Cj and dj {j = 1,2) are defined in Appendix A. 

To proceed, we consider the case k = and make use of the fact that for 7 = 
the stability limit r+ coincides with L2 and remains close for 7 <C 1. We therefore 
expand the value of /i on the stability limit, 

// = /ic + ??^, with 7] = ci7 + C27^ + ..., I7I < 1, (B.4) 

where ~ P^I-'-fi^' — p) is the value of /.j on the curve L2 at fc = 0. The small 
parameter 77 is a function of 7 and must vanish when 7 = 0. Moreover the amplitudes 
of the mixed mode Sj^ on may be expressed as 

Rl = Wll + '>fTl2 ^ , -R2 = ?'20 + 7?^7-22 + ?7^r23 • • • . (B.5) 
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where — \/ A/i/(.s' — p) is the ampUtude of the pure mode on the hne L2. The 
coefficients are determined by substituting the expansions ljB.5|) into H3.5|) and 
solving at successive orders. The number of coefficients r^j that must be kept depends 
on the particular resonance. In fact the expansion (jB.5p must be carried out to 
0(77"^^) for an m : n resonance. In order to calculate the Cj we substitute (jB.SII into 
(|B.2|) and solve order by order. Finally, the wave number Q is obtained by substituting 
these results into (|B.3ll . 

We now sketch the necessary calculations for the case TO:n = 2:5atA: = 0. 
Inserting IIB.4p and (|B.5|) into (|3.5|) yields 



' 11 



2pr2or22 = 1,1^ r\^^ {s' - p)/{ss' - pp') 



(B.6) 



2s'r2or22 + pVf^ = 1. / ^ ^22 = {s - p')/2r2oiss' - pp') 

at 0(?7^), while at 0{r]^) we have 

2sriiri2 + 2pr2or23 = -vr^afii. \ ^ ^12^ -vs'rl^rll},/2{ss' - pp'), 
s'r2or23 + p'riiri2 = J r23 = lyp'rf^r^o^^ /2{ss' - pp'). 

(B.7) 

If n > 5 we would have found ri2 = r23 ~ and it would be necessary to consider 
higher order corrections. Near the curve L2 (i.e., near the point (/i, fc) = (/ic, 0)) (see 
Fig. 14. 1() the functions gj, for j = 0, 1, 2, may be written as 

.90 = ^" (4^V(ss' - pp')^ii^2™o+') + 0(^"+i), 

.91 = (455'(ss' - pp')^?i^2o) + ?7""^(2s'n + A^(n - l)r™)MVar^ (B.8) 

+ -f^d'{-2.s'rl„ - 77'4sV2or22 - 4.s' r2or23) + 0(r/" + 7^""') 
52 = 5'<5'(2sV2o) + 0(?7'), 



where only those terms necessary to compute Q have been kept in (|B.8p . Finally 
composing the expansions ljB.8|l and (jB.4p and substituting the result into (|B.3p we 
find (after some algebra) the result (|4.9() . 
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